Roles of deformation and neutron excess on the giant monopole resonance in 

neutron-rich Zr isotopes 
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We investigate the roles of deformation on the giant monopole resonance (GMR), particularly 
the mixing of the giant quadrupole resonance (GQR) and the effects of the neutron excess in 
the well-deformed nuclei around 110 Zr and in the drip-line nuclei around 140 Zr by means of the 
deformed quasiparticle-random-phase-approximation employing the Skyrme and the local pairing 
energy-density functionals. It is found that the isoscalar (IS) GMR has a two-peak structure, the 
lower peak of which is associated with the mixing between the ISGMR and the K* = + component 
of the ISGQR. The transition strengths of the lower peak of the ISGMR grows as the neutron num- 
ber increases. In the drip-line nuclei, the neutron excitation is dominant over the proton excitation. 
We find for an isovector (IV) excitation the GMR has a four-peak structure due to the mixing of 
the IS and IV modes as well as the mixing of the K n — + component of the IVGQR. Besides the 
GMR, we find the threshold strength generated by neutrons only. 

PACS numbers: 21.10.Re; 21.10.Pc; 21.60.Ev; 21.60.Jz; 24.30.Cz 



I. INTRODUCTION 

Nuclei far from the stability has been one of the major 
interests in nuclear structure physics because exploring 
the multipole responses in unstable nuclei can provide in- 
formation on collective excitations in asymmetric nuclear 
systems. The surface structure of neutron-rich nuclei is 
quite different from that of stable nuclei, because of the 
presence of the loosely bound neutrons. Since collective 
excitations are sensitive to the surface structure, we can 
expect new kinds of exotic excitation modes to appear in 
neutron-rich nuclei. 

An example of such an exotic excitation mode in un- 
stable nuclei is the pygmy dipole resonance (PDR) ap- 
pearing in the low energy region of the giant dipole reso- 
nance (GDR) [J. The PDR consists of the isoscalar (IS) 
and isovector (IV) characters, and they are thought to 
be strongly mixed Q. Strong mixing of the IS and IV 
character leads to a possible emergence of the neutron- 
or proton-dominant giant resonance (GR) in unstable nu- 
clei. 

GRs are typical collective modes of excitation in nu- 
clei. Effects of the nuclear deformation on the GRs has 
been investigated for long both experimentally and the- 
oretically Q. The deformation splitting of the GDR 
has been well established. This is due to the oscillation 
along the major- and minor-axis Q. For the GRs with 
higher multipolarity, the deformation splitting is less pro- 
nounced [3|. Deformation effects on the GMR has been 
discussed soon after the GMR was established [5] . 

Collective excitations such as GRs and low-lying vi- 
brational modes are described microscopically by the 
self-consistent Hartree-Fock (HF) plus the random-phase 
approximation (RPA). There have been many attempts 
to investigate the new kinds of collective excitations in 
unstable nuclei [H, 0|- They are however largely re- 
stricted to spherical systems, and the collective excitation 
modes in deformed neutron-rich nuclei remain mostly un- 



explored. 

Role of deformation on GRs has been studied by means 
of the deformed quasiparticle-RPA (QRPA) employing 
the Gogny interaction in the light mass region [8[ . Using 
Skyrme functionals, GRs in heavy systems have been in- 
vestigated by means of the RPA, where the separable ap- 
proximation is employed for the residual interaction [9j . 

The ground-state properties, deformation and super- 
fluidity, in neutron-rich Zr isotopes up to the drip line 
have been studied QUdlj], and it is predicted that 112 Zr 
has large deformation, and that the drip-line nuclei are 
also strongly deformed. In the present article, we in- 
vestigate the GMR in neutron-rich medium-mass nuclei 
around 110 Zr and the drip-line nuclei around 140 Zr with 
paying attention to the deformation effects and the mix- 
ing of the IS and IV characters. To this end, we use the 
newly-developed parallelized HFB and deformed-QRPA 



calculation code [12[. 



The article is organized as follows: In Sec. HH the de- 
formed Skyrme-HFB + QRPA method is recapitulated. 
In Sec. IIIH results of the numerical analysis of the gi- 
ant resonances in deformed neutron-rich Zr isotopes are 
presented. Finally, summary is given in Sec. IIVI 



II. MODEL 

A detailed discussion of the deformed Skyrme-HFB + 
QRPA can be found in Refs. [l3j]. Therefore, we just 
briefly recall the outline of the calculation scheme. 

In order to describe simultaneously the nuclear defor- 
mation and the pairing correlations including the un- 
bound quasiparticle states, we solve the HFB equa- 
tions jl4fl5j in coordinate space using cylindrical coordi- 
nates r — (p,z,<j)). We assume axial and reflection sym- 
metries. For the mean- field Hamiltonian h, we employ 
the SkM* interaction [l6| . The pairing field is treated by 
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using the density-dependent contact interaction (11 



u pair 



(r,r') = 
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where po( r ) denotes the isoscalar density of the ground 
state and P a the spin exchange operator. Assuming time- 
reversal symmetry and reflection symmetry with respect 
to the x — y plane, we have to solve for positive il and 
positive z only, 51 being the z— component of the angular 
momentum j. We use a lattice mesh size Ap = Az = 0.6 
fm and a box boundary condition at p max — 14.7 fm, 
z max = 14.4 fm. The differential operators are repre- 
sented by use of the 11-point formula of finite difference 
method. Since the parity and f2 are good quantum num- 
bers in the present calculation scheme, we have only to di- 
agonalize the HFB Hamiltonian for each sector. The 
quasiparticle energy is cut off at -Eqp iC ut = 60 MeV and 
the quasiparticle states up to $1* = 23/2 ± are included. 

The pairing strength parameter t' Q is determined so 
as to approximately reproduce the experimental pair- 
ing gap of 120 Sn (A cxp = 1.245 MeV). The strength 
t' Q = —240 MeV fm 3 for the mixed- type interaction 
(t' 3 = -18.75to) 0] with 7 = 1 leads to the neutron 
pairing gap (A v ) = 1.20 MeV in 120 Sn. 

Using the quasiparticle basis obtained as a self- 
consistent solution of the HFB equations, we solve the 
QRPA equation in the matrix formulation 19] 
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The residual interaction in the particle-hole (p-h) channel 
appearing in the QRPA matrices A and B is derived from 
the Skyrme density functional. We neglect the Coulomb 
interaction to reduce the computing time in the QRPA 
calculation. We expect that the two-body Coulomb in- 
teraction plays only a minor role [13,[2l|. We drop the so- 
called "J 2 " term C t T both in the HFB and QRPA calcu- 
lations. The residual interaction in the particle-particle 
(p-p) channel is derived from the pairing functional con- 
structed from the density-dependent contact interaction 
([T]). Furthermore, we introduce a cut-off for the two- 
quasiparticle excitation energy as 60 MeV. 

Since the full self-consistency between the static mean- 
field calculation and the dynamical calculation is broken 
by the above neglected term and the severe cut-off for 
the two-quasiparticle excitations, the spurious states may 
have finite excitation energies. In the present calculation, 
the spurious states appear at 0.58 MeV and 0.53 MeV for 
the — + and 1 + states, respectively in 110 Zr. We 
assume the contamination of the spurious component in 
GRs to be small because the GRs are well apart from the 
spurious states. 

In the present calculation, the number of two- 
quasiparticle excitations for the K v = + excitation in 
no Zr (i40 Zr ) is about 36 ( 46 200). For the calcula- 
tion of the matrix elements and diagonalization, it takes 



TABLE I: Ground state properties of 
100,108,110,112.136,138,140 Zr obtained by the de f rmed HFB 

calculation with the SkM* interaction and the mixed-type 
pairing interaction. Chemical potentials (in MeV), defor- 
mation parameters, intrinsic quadrupole moments (in fm 2 ), 
average pairing gaps (in MeV), root-mean-square radii (in 
fm) for neutrons and protons are listed. 





100 Zr 


108 & 


110 Zr 


112 Zr 


136 Zr 


138 & 


i4o Zr 


\ v 


-6.81 


-4.57 


-4.31 


-3.76 


-0.51 


-0.36 


-0.10 


A-tt 


-10.6 


-13.6 


-14.4 


-15.0 


-22.3 


-22.6 


-22.8 


P2 


0.38 


0.37 


0.37 


0.39 


0.28 


0.29 


0.31 


02 


0.41 


0.42 


0.42 


0.43 


0.32 


0.33 


0.36 


(Q)u 


602 


791 
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947 


1023 


1135 


(Qh 


403 


429 


432 


456 


355 


374 


409 


(A)* 


0.00 


0.61 


0.55 


0.00 


0.70 


0.64 


0.53 


(A). 


0.32 


0.00 


0.00 


0.00 


0.36 


0.36 


0.36 




4.60 


4.77 


4.82 


4.88 


5.29 


5.34 


5.39 




4.43 


4.51 


4.53 


4.56 


4.69 


4.72 


4.75 



about 346 (570) CPU hours and 184 (310) CPU hours, 
respectively. 



III. RESULTS AND DISCUSSION 

We summarize in Table Q] the ground state properties. 
The neutron-rich Zr isotopes under investigation are pro- 
lately deformed. This is consistent with the results calcu- 
lated using the Skyrme SLy4 interaction [l(| , and we also 
find that 112 Zr has largest deformation. Both neutrons 
and protons are unpaired in this nucleus. This calcula- 
tion indicates that 140 Zr is located very close to the drip 
line; the chemical potential is only —0.1 MeV. The neu- 
tron skin well develops in this nucleus; the difference of 
the radii of neutrons and protons is 0.64 fm. 

Figures [T] and [2] show the isoscalar (IS) and isovector 
(IV) monopole and quadrupole transition-strength dis- 
tributions in 90,100,108,110,112 Zr and m 136,138,140 Zr Here 

we define these operators as 



1-2=1,-1 j 

=1,-1 "* 

Y, I drr 2 Y XK (f)i;l(r)^(r), 
1, =1,-1 

H=Ik = E / drr z r 2 Y XK (f (r)j> Tz (r) 



pr=l 



PT=0 



(3) 
(4) 
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A nucleon creation operator ip\r) at the position r is 
written in terms of the quasiparticle wavefunctions. 
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FIG. 1: (Color online) Isoscalar (IS) and isovector (IV) monopole and quadrupole transition strengths in smoo, 108,110,112 Zr. The 
neutron and proton transition strengths are also shown for the monopole excitation. The arrows indicate the one-quasineutron 
continuum threshold E t h,in = |A| + min.Ec, and the two-quasineutron continuum threshold E t h,2n = 2|A|. Shown is only the 
one-neutron continuum threshold for 90 > 100 > 112 Zr, the normalfluid systems. 



The strength function is calculated as 



= EE 



r/2 



(i\F- 
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The smearing width T is used to remove the spurious os- 
cillation in the strength functions due to the box bound- 
ary condition, and it is set to T = 2 MeV in the case of 
the box size of 15 fm [22j |. 

The ISGMR has a two-peak structure in the deformed 
nuclei. The higher-energy peak of the IS monopole ex- 
citation seen around 20 MeV is identified as a primal 
ISGMR because the ISGMR in 90 Zr is located at around 
18 MeV, where we have only one peak. 

To understand the origin of the lower-energy peak 
around 10 MeV, we show in Figs. [DM-EHo) the K n = 0+ 
component of the quadrupole transition strength. At 
the same energy region where the lower-energy peak 
of the ISGMR appears, we can see a peak structure. 
It is noted here that the neutron-rich Zr isotopes un- 
der investigation are prolately deformed, so the giant 
quadrupole resonance (GQR) splits into three resonances 
of the = + ,l + , and 2 + components, and the reso- 
nance peak of the K T = + component lowers in energy. 

We are going to discuss in detail the deformation split- 
ting of the ISGMR. Figure EJa) shows the mean energy 
of the lower peak and the upper peak of the ISGMR. 



The mean energy is evaluated as mi /mo, the ratio of the 
energy-weighted sum to the non energy-weighted sum. 
The energy interval is set to 5 MeV< haji < E c and 
E c < hcji < 25 MeV for the lower and upper peaks. The 
value E c is set to 15 MeV and 14 MeV for 100 Zr and 
108,110:112 Zr, respectively. We change E c associated with 
the mass number because the resonance energy decreases 
as the mass number increases. The values obtained in the 
scaling model [23| are also shown in Fig. [3] 

The results of the present calculation are not very dif- 
ferent from those obtained in the scaling model, with the 
exception that the scaling model overestimates the en- 
ergy of the compressible mode [23]. As the system is 
deformed from 90 Zr to 100 Zr, the excitation energy of the 
upper peak is upward-shifted due to the coupling with 
the K n = 0+ component of the ISGQR. Then, as the 
mass number increases from A = 100 to 112, the exci- 
tation energy is downward-shifted. The difference in the 
excitation energy of the upper peak in 100 Zr and in 112 Zr 
is 0.7 MeV. This energy difference is understood by the 
systematic trend of the excitation energy E x oc A -1 / 3 . 

The lower peak of the ISGMR in deformed systems 
originates from the ISGQR in the spherical limit. The 
excitation energy of the lower peak in the scaling model 
agrees well with that obtained in the present calculation. 

Figure [SJb) shows the fraction of the EWSR value 
for the lower and upper peaks of the ISGMR. The en- 
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FIG. 2: (Color online) Same as Fig. [T] but in 
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FIG. 3: (Color online) (a) Mean energy of the ISGMR in 
Zr isotopes. The energy of the upper and lower peaks are 
denoted by triangle and square, respectively, (b) Fraction of 
the EWSR value for the upper and lower peaks. The mean 
energy of the upper and lower peaks is calculated as mi /mo 
for the energy interval of 5 MeV < E < E c and E c < E < 25 
MeV, and the fraction is calculated for each interval. The 
lines represent the results obtained in the scaling model [23l ]. 



ergy interval to evaluate the energy weighted strength 
is the same as in the evaluation of the mean excitation 



energy. The values obtained in the scaling model 23 1 
are also shown. Unlike for the excitation energy, we can 
clearly see deviation from the values obtained in the scal- 
ing model as the neutron number increases. We have an 
appreciable enhancement of the transition strengths of 
the lower peak. One may attribute the deviation to the 
neutron excess. 

To clearly see the effect of the neutron excess, we move 
on to the drip-line systems. Figure [2] shows the IS and 
IV monopole and quadrupole transition-strength distri- 
butions in 136 >i38,i402 r> can see t ne further enhance- 
ment of the IS monopole transition strength in the lower 
peak region. The mean energy of the lower and upper 
peaks is shown in Fig. (3Ja). We set E c to 13 MeV. Even 
in the drip-line nuclei, we can not see an appreciable de- 
viation from the simple scaling model both for the mean 
energies of the lower and upper peaks. 

However, the fraction of the EWSR value of the lower 
peak keeps increasing toward a drip line. Note that 
the deformation of near-drip-line nuclei 136 . 138 > 140 Zr is 
smaller than that of 108 ' 110 . 112 z r . In the scaling model, 
the coupling between the GMR and GQR is governed 
by the quadrupole deformation of the ground state [23| . 
Thus, the calculated fraction of the EWSR value of the 
lower peak in the scaling model is small in the A ~ 140 
region. 

To investigate the origin of the enhancement of the 



TABLE II: Sum of the transition strengths of neutrons and 
protons Ml, Ml (in 10 3 fm 4 ) and M„/M^ divided by the 
neutron and proton numbers for the lower and upper peaks 
of the ISGMR. 



Lower peak 



Upper peak 



Ml Ml 


M v /M w 
N/Z 


Ml 


Ml 


M u /M* 
N/Z 


90 'Zr 




2.62 


1.70 


0.99 


100 Zr 2.33 0.75 


1.18 


2.75 


1.44 


0.92 


108 Zr 3.53 0.76 


1.27 


3.68 


1.43 


0.94 


110 Zr 3.64 0.77 


1.24 


3.85 


1.48 


0.92 


112 Zr 4.50 0.85 


1.28 


3.93 


1.54 


0.89 


136 Zr 9.35 0.60 


1.64 


5.82 


1.78 


0.75 


138 Zr 10.2 0.65 


1.62 


5.83 


1.82 


0.73 


140 Zr 10.8 0.72 


1.55 


6.05 


1.86 


0.72 



transition strengths in the low energy region, we show 
in the middle panels of Figs. [T] and [2] the neutron and 
proton monopole transition strengths. It is seen that the 
neutron excitation plays a dominant role in increasing 
the IS transition strengths in the energy region around 
10 MeV. 

Table [TT] summarizes the contribution of the neutrons 
and protons to the lower and upper peaks of the ISGMR. 
We can see the effect of the neutron excess on the GMR 
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quantitatively from this Table. Listed in this Table are 
the sum of the transition strengths of neutrons M 2 . and 
protons M 2 . in the energy intervals of 5 MeV < fiuJi < E c 
and E c < hwi < 25 MeV. If the neutrons and protons 
coherently contribute to the excitation mode, the ratio of 
the matrix elements of neutrons to protons may approach 
the ratio of the neutron number to the proton number. 

The ratio of the matrix elements divided by the ra- 
tio of neutron to proton numbers, {M u /M 7r )/{N/Z), for 
the lower peak increases in the drip line nuclei. On the 
other hand, the quantity (M JJ /M 7r )/(N/Z) for the upper 
peak decreases in the drip line nuclei. Thus, the neutron 
excitation concentrates particularly in the energy region 
of the lower peak of the ISGMR. Besides the GMR, we 
can see appearance of the threshold strength consisting 
of neutron excitations only in the drip- line nuclei. This 
fact is similar to the findings of Ref. [24] in 60 Ca. 

We are going to move on to the discussion on the prop- 
erties of the IVGMR. In 90 Zr, the mean energy of the 
IVGMR is 31.6 MeV. The hydrodynamical model taking 
a surface effect into account [25| predicts the excitation 
energy of the IVGMR to be 

S a ~88A- 1 / 6 (l + — A" 1 / 3 )- 1 / 2 MeV. (8) 
3 

This gives the energy of the IVGMR in 90 Zr to be about 
29 MeV. This is consistent with the microscopically cal- 
culated energy. 

The IVGMR in deformed neutron-rich nuclei has inter- 
esting features, where a four-peak structure appears as 
shown in Fig. [TJ The four-peak structure is clearly seen 
especially in 112 Zr. 

The highest peak above 30 MeV corresponds to a pri- 
mal IVGMR by comparing the peak energy with that in 
90 Zr. We have a peak structure at 22 — 23 MeV where the 
IVGQR is seen. Thus the second highest peak emerges 
associated with the coupling to the — + component 
of the IVGQR. 

Besides these two peaks, we can see two more peaks 
below 20 MeV. We can consider that they are brought to 
the energy region of the ISGMR due to the neutron ex- 
cess. One is associated with the coupling to the ISGMR, 
and the other is due to the coupling to the K* = + 
component of the ISGQR. 

In the drip-line nuclei, the transition strengths concen- 
trate in the ISGMR region. We can still see a four-peak 
structure in 13 6,i38,i402 r> ^he jy transition strengths be- 
come enhanced in the energy region where the neutron 
transition strengths increase. Therefore, we can see the 
threshold strength in the IV monopole strength distribu- 
tion as well. 



IV. SUMMARY 

We have investigated the deformation effects on the 
GMR by employing the newly developed parallelized 
QRPA calculation scheme. Associated with the defor- 
mation, the excitation modes with the angular momenta 
I = and / = 2 can mix in the K v = + channel, and 
accordingly we have obtained a two-peak structure of the 
ISGMR. The lower peak of the ISGMR is due to the cou- 
pling with the ISGQR. The deformation splitting is seen 
in nuclei from near-stability line to drip line. 

The deformation splitting obtained in the present cal- 
culation is in fairly good agreement with the result of the 
scaling model, while the distribution of the strengths is 
quite different when the drip line is approached. The IS 
monopole transition strengths are enhanced in the energy 
region of the lower peak of the ISGMR. This is due to 
the concentration of the neutron transition strengths in 
the lower energy region. 

We have found a complicated peak structure for the 
IVGMR in deformed neutron-rich nuclei, where a four- 
peak structure appears. The highest peak corresponds 
to a primal IVGMR, and the second is associated with 
the coupling to the — + component of the IVGQR, 
and the two more peaks are found in the ISGMR region 
associated with the neutron excess. 

Besides the GMRs, the threshold strength is found in 
the drip- line nuclei. The strength consists of neutron 
strengths only. Since the low-lying IC = + mode is 
quite sensitive to the shell structure and surface struc- 
ture [26( | . we are going to investigate the microscopic 
structure of the low-lying X* = + states with special 
attention to the coupling to the continuum states and the 
pairing correlations. 
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